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Abstract 

We generate inherent structures, local potential-energy minima, of the "/c-space overlap poten- 
tial" in two-dimensional many-particle systems using a cooling and quenching simulation technique. 
The ground states associated with the /c-space overlap potential are stealthy (i.e., completely sup- 
press single scattering of radiation for a range of wavelengths) and hyperuniform (i.e., infinite 
wavelength density fluctuations vanish). However, we show via quantitative metrics that the in- 
herent structures exhibit a range of stealthiness and hyper uniformity depending on the fraction 
of degrees of freedom that are constrained. Inherent structures in two dimensions typically con- 
tain five-particle rings, wavy grain boundaries, and vacancy-interstitial defects. The structural 
and thermodynamic properties of inherent structures are relatively insensitive to the temperature 
from which they are sampled, signifying that the energy landscape is relatively flat and devoid 
of deep wells. Using the nudged-elastic-band algorithm, we construct paths from ground-state 
configurations to inherent structures and identify the transition points between them. In addition, 
we use point patterns generated from a random sequential addition (RSA) of hard disks, which 
are nearly stealthy, and examine the particle rearrangements necessary to make the configurations 
absolutely stealthy. We introduce a configurational proximity metric to show that only small local, 
but collective, particle rearrangements are needed to drive initial RSA configurations to stealthy 
disordered ground states. These results lead to a more complete understanding of the unusual be- 
haviors exhibited by the family of "collective-coordinate" potentials to which the fc-space overlap 
potential belongs. 
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I. INTRODUCTION 



Recently, we have been interested in a family of soft, long-range pair interactions that 
give rise to novel physical behaviors of many-particle systems including classical disordered 
ground states for a range of densities, negative thermal expansion, and vanishing normal- 
mode frequencies.- - -^- This family of soft pair interactions includes those pair potentials 
v(r) with Fourier transforms V(k) that are positive, bounded, and vanish to zero at some 
finite wavenumber K and beyond. The potential energy $ of a system of N particles in a 
fundamental cell with volume Q under periodic boundary conditions can be written as 



where = [r» — r 3 -| is the distance between particles i and j, p(k) are the Fourier coeffi- 
cients of the density field, and k are the wave vectors appropriate for the system size and 
shape. These soft pair interactions, called "collective-coordinate" potentials, are bounded 
and possess long-range oscillations in real space.- 

Because of the finite cutoff in /c-space, an analytic lower bound for the potential energy 
per particle is easily obtainable, and numerical methods allow for the construction and 
investigation of ground states with extremely high precision. A series of studies has examined 
the structural characteristics of ground states in one,- two,- and three dimensions.- Three 
structural ground-state regimes were found to characterize these two-dimensional systems. 
Increasing the fraction of degrees of freedom that are constrained \i equivalently decreasing 
density, spans regimes that are disordered, wavy- crystalline, and crystalline. Recently, an 
analytical connection between the fraction of constrained degrees of freedom within the 
system and the disorder-order phase transition for a class of target structure factors has 
been provided 9 by examining the realizability of the constrained contribution to the pair 
correlation function. 

More recently, these potentials were used to construct "stealthy" and "hyper uniform" 
materials in two dimensions.- "Stealthy" materials refer to point patterns that completely 
suppress single scattering of radiation for certain wavelength ranges. While nonstealthy 
materials may allow for some small amount of scattering, stealthy materials are absolutely 
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FIG. 1: (Color online) Schematics of the equilibrium equation of state (solid lines) and the nonequi- 
librium path that results in glassy behavior (dotted lines) along an isobar. With familiar potentials 
(e.g. Lennard- Jones) (left), the equilibrium curve shows a first-order phase transition while the 
glassy curve undergoes supercooling and a glass transition. With the soft, long range fc-space 
overlap potential (right), high-pressure systems result in a first-order transition that is not present 
for low-pressure systems. Both collective-coordinate systems undergo negative thermal expansion. 

transparent at those wavelengths. "Hyperuniform" refers to point patterns in which infinite 
wavelength density fluctuations vanish. 10 (These terms are defined more precisely in Sec. 
UT1 ) The hyper uniformity notion enables the rank ordering of crystals, quasicrystals and 
special disorderedmany-particle systems.— ~— Interestingly, disordered hyperuniform many- 
particle ground states, and therefore also point distributions, with substantial clustering can 
be constructed.- 

Constructed ground states are now the basis for novel disordered materials with tunable, 
photonic band gaps.— ^ Using a model potential from the family of collective-coordinate po- 
tentials in two dimensions, the "fc-space overlap potential," described in Sec. [TTJ we observed 
negative thermal expansion and vanishing normal-mode frequencies for ground-state config- 
urations. We attributed these phenomena to the nature of the underlying energy landscape, 
which we described as having ground-state "valleys" that weave between the higher-energy 
portions of the landscape.- 1 ^ 

Despite the discovery of the existence of such unique properties, the fundamental mecha- 
nisms allowing for disordered ground states and stealthy point patterns still need to be fully 
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elucidated. Here, we analyze the energy landscape further by examining the associated in- 
herent structures, or local potential-energy minima, for the aforementioned /c-space overlap 
potential in two dimensions to better understand these properties. 

We have two primary motivations for this work. First, through the course of our previous 
research, we have observed several paradoxical phenomena related to local potential-energy 
minima and glassy behavior. While an "inherent structure" is in general any local potential- 
energy minimum, the term "glass" refers to an amorphous solid that is kinetically trapped in 
a potential-energy well. In Ref. [2], we observed an unusual equation of state when cooling 
a system compared to a more familiar (e.g., Lennard- Jones) glass-forming system as Fig. 



[T] demonstrates this schematically. In Fig. 1(a), the dark and dotted lines represent the 



equilibrium equation of state and nonequilibrium isobaric cooling path for a Lennard- Jones 



system. Figure 1(b) shows the corresponding cooling paths for a system interacting with 
the fc-space overlap potential. In the more familiar glass-forming system, there is a well- 
defined glass transition below the freezing point that is dependent on the rate at which 
the system is cooled. 15 It is well known that the glass- forming behavior of a system is a 
function of the underlying energy landscape.— The depths of the potential-energy wells in 
the energy landscape dictate the dynamics of structural rearrangements of the system.— 1 ^ 
For the fc-space overlap potential, we observed that at high pressure the nonequilibrium 
curve deviated from the equilibrium curve at a temperature above the melting temperature, 



as demonstrated in Fig. 1(b) At lower pressures where there is no well-defined melting 
temperature, as in the wavy-crystalline regime, the nonequilibrium curve deviates similarly 
from the equilibrium curve. In addition, the fc-space overlap potential gives rise to negative 
thermal expansion, generally for nondimensional T* < 0.0007.-^ 

Other paradoxical behaviors have led us toward a broad analysis of the energy landscape. 
We have observed that some ground-state configurations, particularly crystalline structures, 
had many zero-energy normal modes. These ground states are therefore not mechanically 
rigid. However, while using numerical minimization algorithms to construct ground-state 
configurations, we encountered many mechanically stable local potential-energy minima. 
Lastly, there exists a large density range in which the energy landscape was evidently devoid 
of all local minima. An increment in x just outside this range introduced local potential- 
energy minima to the landscape. 

While we have previously observed these phenomena,- - - the fundamental mechanisms 
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that underlie them are not fully understood. The inherent-structure analysis has proven 
to be a fruitful method for relating the energy landscape to low-temperature phenomena. 
Many studies have examined inherent-structure characteristics in glass-forming liquids with 
strong repulsive cores such as the Stillinger- Weber potential,— r — water-like pair potentials,— 
binary Lennard- Jones-like systems,—"— and general repulsive potentials.— 1 ^ The inherent 
structures for Lennard- Jones and steeply repulsive potentials are in general not hyperuniform 
or stealthy due to the dominance of grain boundaries and vacancy defects. 

In contrast to the above list of strongly repelling potentials, the /c-space overlap potential 
is a soft interaction. Soft interactions are often useful models for soft-matter systems such 
as colloids, polymers, and microemulsions. 2 - In addition, these fc-space overlap interactions 
are also qualitatively similar to Friedel oscillations in molten metals.— The fc-space over- 
lap potential is localized in fc-space and delocalized in real space as a result of the Fourier 
transform. Certain duality relations link the ground-state energies of the /c-space overlap 
potential to the ground-state energies of real space analog.- 1 ^ There have been several in- 
vestigations of the inherent structures of various soft interactions. The Gaussian-core model 
in two dimensions has polycrystalline inherent structure with a large correlation length even 
when sampled from the liquid state.™ Energy landscape analyses revealed that the range of 
the Morse potential affects the relation between temperature and the potential-energy distri- 
bution of sampled inherent structures.— While for the Yukawa potential, inherent structures 
varied depending on whether they were obtained from the liquid, crystal or hexatic phase.— 

Novel, stealthy dielectric materials are currently being fabricated,™ and nearly stealthy 
ceramic materials are of interest for optical applications.— 1 ^ Because of these recent ex- 
perimental applications and the unusual physical behaviors discovered, -~- we are also mo- 
tivated to understand the fundamental differences between point patterns that absolutely 
suppress scattering for certain wavelengths and those that nearly suppress scattering. Our 
ground-state construction procedure^ automatically distinguishes stealthy and hyperuniform 
configurations from those that do not possess such properties. In particular, we want to un- 
derstand whether the particle rearrangements required to transform nearly hyperuniform 
and stealthy materials to configurations that are perfectly hyperuniform and stealthy are 
global or local in nature. 

While a general method for understanding these particle rearrangements would require 
a new algorithm to search for collective motions that increase stealthiness, the collective- 
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coordinate approach provides an excellent framework from which to address this question. 
Here, we provide two methods for identifying rearrangements in particle systems to achieve 
stealthy and hyperuniform materials. For large Xi we use the nudged-elastic-band algorithm 
to connect inherent structures, which are nearly hyperuniform and nearly stealthy, to ground 
states along a minimum-energy path. For small x values, we have identified that there 
are no inherent structures higher up in the energy landscape. Therefore, we study the 
rearrangements from a saturated random sequential addition (RSA) of hard disks 37 to ground 
states of the fc-space overlap potential. We use saturated RSA systems as initial conditions 
because upon saturation they suppress scattering for small wavenumbers and are nearly 
hyperuniform.— - — We also introduce a stealthiness metric and a configurational proximity 
metric to quantify characteristics of these transitions. 

In this paper, we use a collective-coordinate potential and a simulation methodology to 
find stealthy point patterns and local potential-energy minima above the ground state. We 
probe the following fundamental questions: 

• How are inherent structures and their thermodynamic properties in collective coordi- 
nate systems different from those found in other soft-potential systems? 

• To what extent are inherent structures stealthy and/or hyperuniform? 

• How are the features of inherent structures related to the pair potential function? 

• What collective particle rearrangements are necessary to construct a path from an 
inherent structure to a ground state? 

• What global and/or local particle rearrangements are necessary to convert a nearly 
stealthy and nearly hyperuniform system in to one that is absolutely stealthy and 
hyperuniform? 

The remainder of this paper is organized as follows. In Sec. [Til we define the fc-space over- 
lap potential, introduce several definitions, and briefly review the ground-state structural 
regimes. In Sec. 1111} we detail the methods we use to obtain and characterize inherent struc- 
tures. In addition, we discuss the methods we use to find transition states between inherent 
structures and ground states. In Sec. IIV| we examine the thermodynamic properties of in- 
herent structures as a function of system size and cooling rate. The inherent structures are 
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characterized in Sec. |VJ Paths connecting inherent structures to ground states are provided 
in Sec. I VI I while the rearrangements from RSA systems to ground states are explored in Sec. 
IVHI Concluding remarks are provided in Sec. IVIIII 



II. COLLECTIVE COORDINATES AND THE OVERLAP POTENTIAL 

While a fully detailed summary of collective coordinates is provided in Ref . [2] , we provide 
a brief summary of the relevant mathematical relations and definitions here. Recall that the 
potential energy of the system of interacting particles in a periodic simulation box is defined 
in Eq. ([T]). The collective density variables p(k) are given by 

N 
3=1 

where rj is the position of the j th particle. The wave vectors for a periodic box correspond 
to linear combinations of the reciprocal lattice vectors associated with the periodic box. For 
a rectangular box of dimensions L x and L y , those wave vectors are k = [2irn x /L x , 27rn y /L y ], 
where n x and n y are integers. The structure factor 5*(k), proportional to the intensity of 
scattering of radiation, is related to the collective density variables via 

W - <«. (4) 

where < • • • > is the appropriate ensemble average. When angularly averaged, the structure 
factor S(k) is only dependent on the wavenumber k = |k|. Since the collective-coordinate 
class of potentials are bounded positive for k < K and zero for k > K, any configuration 
for which S(k) is constrained to be zero for all < |k| < K is a ground state. 
For such potentials, it is useful to refer to the dimensionless parameter 

M(K) 

as the fraction of degrees of freedom that are constrained, where M(K) is the number of 
independent wave vectors < |k| < K and dN is the total degrees of freedom.- In two 
dimensions \ is limited to the range \ < 0.91. This limit was discussed in Ref. [2] and 
arises from the inability to suppress Bragg scattering of the triangular lattice.— Here, we 
fix K = 1, and therefore x is inversely related to the number density p = N/Q. Some 
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FIG. 2: The overlap pair potential function v(r) in the infinite-volume limit and its corresponding 
Fourier transform V(k), the /c-space overlap potential, as adapted from Ref. [2]. 

relations between \ and p are given in Ref. [2]. Three structural regimes exist for ground- 
state configurations in two dimensions. For x < 0.577, ground states are disordered. For 
0.577 < x < 0.780, they are wavy crystalline, while for 0.780 < x < 0.91, crystalline 
structures are the most viable ground-states configurations identified via numerical search 
techniques.-^ We have observed that quenching systems with x < 0.5 always resulted in 
a ground state. The energy landscape is evidently devoid of local potential-energy minima 
above the ground state when x < 0.5 and thus inherent-structure analysis is limited to those 
systems with x > 0.5. 

We continue to use the fc-space overlap potential, introduced by Torquato and Stillinger, 5 
since it is exactly representable in real and reciprocal space, and is sufficiently short-ranged 
in real space for computational purposes; see Fig. |2j The fc-space overlap potential is 
proportional to the intersection area between two disks of diameter K with centers separated 
by k, and hence for k < K is 

k 2 ^ 1/2 



V(k) = ^ 

7T 



COS 



k \ k 

KJ~KV K 2 



(6) 



and zero for all k > K. In the infinite- volume limit, the associated real space pair potential 
is 

\L T / Kr \ 1 2 

(7) 



( \ v ° 



where J± is the Bessel function. Henceforth, we refer to this v(r) as the "overlap potential." 
The potential v(r) is bounded at r = and behaves as cos 2 (Kr/ '2 — 37r/4)/r 3 for large r. 



In this paper, we choose to report a scaled potential energy e, omitting the structure- 
independent contribution at k = and removing the additive constants, i.e., 

so that ground state energies are equal to zero regardless of x f° r X < 0.91. This new scaled 
energy provides a sense as to how much a configuration differs energetically from a ground 
state, and is related to the actual potential energy, choosing V$ = 1, via 

</> = $/iv = £ WMCQ*)] ( 9 ) 

k 

= f"^Enk) + e (10) 

k 

Note that the structure-independent terms sum to be nonnegative because v(r) is a nonneg- 
ative function. We choose the cutoff for ground states to be those configurations in which 
e < 10- 10 . 

We also compute a stealthiness metric r\ as the average of the values of the structure 
factor for all |k| < K, defined to be 

V = S(k)/M(K) (11) 

o<|k|<_ft: 

This metric is a measure of the stealthiness because it vanishes to zero for completely 
transparent systems but will remain nonzero for nonstealthy systems. While there is a 
numerical limitation of 10 -16 for double precision computing, there are certain practical 
limits on stealthiness. Because of the weights assigned by V(k) for k near K, the metric rj 
may not necessarily be suppressed to zero for ground states due to numerical precision of the 
quantity V(k)S(k). Therefore, the stealthiness metric is best used as a relative comparison 
of stealthiness between configurations and not as a determination of a configuration as a 
ground state. 

III. METHODS 

Our method to generate inherent structures follows that of Sastry and coworkers,— which 
is based on the initial algorithms of Stillinger and Weber.— Initially, we simulate the 
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atomic motions of the system interacting via the fc-space overlap potential in Eq. (GQ) by 
integrating the Newtonian equations of motions, assuming unit mass, 

d 2 Ti 



-V<f>. (12) 



dt 2 

We use the velocity Verlet algorithm^ with a time step 5t of 0.4, which is chosen so as to 
accurately conserve total energy in the constant NVE ensemble. Throughout the dynamical 
simulation, configurations were sampled and subjected to a quenching of the potential energy. 
Here, the sampled configuration is used as the initial condition for a conjugate gradient 
minimization. 

Upon termination of the conjugate gradient quenching, we subject the system to an addi- 
tional quenching via the MINOP algorithm.— & We previously reported the efficiency of the 
MINOP algorithm compared to the conjugate gradient method.- 1 ^ We add this final quench 
for reasons of increased efficiency and precision compared to simply using the conjugate 
gradient method to high precision. Many of the inherent structures that we have found have 
very slight differences in the scaled potential energy e that can be on the order of magnitude 
10~ 8 . Our implementation of the conjugate gradient method becomes inefficient at these 
tolerances. While there is a possibility that the MINOP algorithm can traverse to other cap- 
ture basins, we suspect that the system is quenched sufficiently deep into a potential-energy 
well after the conjugate gradient step that the MINOP algorithm will not push the system 
into another capture basin. 

We initialize the dynamical trajectories at a temperature associated with the liquid state. 
Here, we report the dimensionless temperature T* = fe^T/Vo, where k% is the Boltzmann 
constant and T is the simulation temperature. We previously reported that the transitions 
from ground states to highly disordered liquids occur in the range of T* = 0.0003-0.001, 
depending on x-~ We initialize our simulations at T* > 0.0018. Configurations are sampled 
every n s time steps until at least fifty are obtained. Then the velocities of the system are 
rescaled to a lower temperature by dT*, typically a value of 0.00015, and the configurations 
are sampled again at the reduced temperature. The rate at which the system is effectively 
cooled becomes 7 = dT* / (50n s St) . We have explored a range of coolings rates from 10 -6 < 
7 < 10~ 9 with n s varying from 3 to 750. Surprisingly, as discussed in Sec. IIVI while the 
properties of the thermal structures are sensitive to the cooling rate, the properties of the 
inherent structures are not. 
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To construct pathways between ground-state configurations and inherent structures for 
X > 0.5, we use the "nudged-elastic-band algorithm" used for finding minimum-energy tran- 
sition paths.— & This method requires as an input two local minima that are nearby in the 
energy landscape. We generate this pair of structures by initializing a known ground-state 
configuration (obtained by the search methods in Ref. [3]), assigning a small temperature 
to the system, and integrating the equations of motion. The configurations are quenched 
by the conjugate gradient/MINOP minimization every five time steps. A simulation run 
is terminated when the quenching of a configuration produces a local minimum in a new 
capture basin. This new inherent structure is presumed to be near the original ground state 
in configuration space. 

We then apply the nudged-elastic-band (NEB) algorithm^^ to the ground state/inherent 
structure pair. The algorithm discretizes the path between two local minima in a potential- 
energy landscape into "beads," or image configurations. These beads, representing config- 
urations, are connected via harmonic springs of zero natural length. Then a minimization 
algorithm is applied to minimize simultaneously the force parallel to the string of beads 
and the force perpendicular to the true force. The algorithm bends the string of beads 
around the "hills" in the energy landscape. The end result is a discrete path that provides 
the minimum-energy path from one potential-energy minimum, through a low-order saddle 
point, and to the other potential-energy minimum. While a large number of "beads" allows 
for a more refined minimum-energy pathway, we find that using fifty beads serves our pur- 
poses well by yielding a smooth transition pathway and a reasonable approximation to the 
saddle point. 

Occasionally, there are other inherent structures in the vicinity of the straight-line inter- 
polation between the inherent structure and ground state used as an initial condition. In 
these cases, the NEB algorithm will not produce a smooth transition from the ground state 
to the original inherent structure, and therefore we refine our initial ground state/inherent 
structure pair. For example, when using the triangular lattice as an initial condition for 
say x — 0.6004, there are number of zero-energy modes that the system can traverse during 
the molecular dynamics trajectories. When applying the NEB algorithm to a triangular lat- 
tice/inherent structure pair, the pathway has a preference to bend toward the zero-energy 
path before moving uphill toward the saddle point. This caused the energy of the transition- 
state pathway to appear rugged since the string of beads was pushed near a different inherent 
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structure. In these cases, we searched the string of beads for a new ground-state/inherent 
structure pair. We minimized the potential energy of each bead until we found a new ground 
state/inherent structure pair closer together than the original pair and applied the NEB al- 
gorithm again. This method does not guarantee that we find the closest inherent structure 
to a ground state nor does it guarantee that the saddle point is a first-order (single negative 
eigenvalue) saddle. 

For x < 0.5, we examine particle rearrangements from RSA configurations to ground 
states. To generate RSA packings, disks of diameter D are randomly, sequentially, and irre- 
versibly placed inside a simulation cell so that they do not overlap with any other particles. 
The process is saturated when additional particles cannot be placed without overlapping 
other particles. In practice, this occurs after a very large number of attempted particle 
insertions are rejected. Here, we terminate the process after 1.5 million attempted particle 
placements. Since K is fixed to be unity, we must select a target number density p t that 
corresponds to our target x- Using a desired number of particles in a saturated system of 
750 particles, we then identify the appropriate system area. The simulation cell is chosen to 
be a square box. We can then assign the particles a diameter D so that upon saturation of 
the RSA process, we achieve close to the desired number density. For RSA processes in two 
dimensions, the saturation packing fraction f s is 0.547.— The assigned diameter becomes 
D = 2{f s / , Kp t ) 1 / 2 . In our RSA patters, the number of particles did not match the target 
number of particles due to finite system effects. Given the dimensions of the box and the 
number of particles, we can then assign the appropriate wave vectors and x value to the 
system. The potential energy of the system is then quenched using the conjugate gradient 
method followed by a quenching using the MINOP algorithm as was detailed above. 

To characterize inherent structures, we use several quantitative metrics. The structure 
factor S(k), cf. Eq. provides structural information on the long-ranged ordering of 
the system. Defined above, e is a scaled potential energy and rj is a metric for the extent 
of stealthiness. In addition, we use the small- features of S(k) to compare the extent 
of hyperuniformity between structures. Hyperuniform systems have the feature that S(k) 
vanishes as k approaches zero. We fit the small- regime (i.e. k < 0.5) to a log-polynomial 
equation 

log S(k) = C + dk + C 2 k 2 + C 3 k 3 . (13) 
The parameter Cq can be considered to be a metric for the extent of hyperuniformity, where a 

13 



value that diverges to — oo indicates that a system is hyperuniform. This metric is best used 
for relative comparisons among systems since numerical precision is limited. For example, 
one could say system A is more hyperuniform than system B because it has a lower value 
of Cq. However, Cq is not be a sufficient test to determine if a configuration is absolutely 
hyperuniform (i.e., one should refrain from saying system A is absolutely hyperuniform). 
We also employ the bond-order parameter defined as 



where ^ is the angle between two particles with respect to a fixed, but arbitrary, coordinate 
axis, to quantify the local orientation order. Particles are considered "bonded" if < 10, 
chosen so as to include nearest-neighbors. For the perfect triangular lattice, has a value 
of unity while for the ideal gas \&6 vanishes. Wavy crystals have value of that can range 
from about 0.4 to 0.8. 

We introduce a configurational proximity parameter p defined as 



where r o i and r/^ are the positions of the RSA configuration and ground state configuration 
respectively and r^jy is the mean- nearest-neighbor in the RSA configuration. This metric 
p quantifies the closeness of two different particle configurations. In what follows, we will 
apply it to characterize the relative particle displacements required for a configuration to 
relax to the ground-state energy. 

IV. RESULTS: COOLING AND QUENCHING 

Before discussing the results of the quenching and cooling simulations, we briefly discuss 
some aspects of system-size effects related to the inherent structures of this system. While 
we previously showed that system-size effects were negligible when characterizing the ground 
states^ and the equilibrium properties, 2 we find that there are appreciable system-size effects 
when considering the inherent structures and energy landscape. 

We have compared inherent structures consisting 418 particles to 2340 particles sampled 
from the same temperature. In general, larger systems provide a slightly deeper quenching 
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of the potential energy. In addition, the bond-order parameters for larger systems tend 
to be less than that of smaller systems. This is expected since larger systems often have 
various polycrystalline-like domains that impart destructive interference on \l/ 6 . The fraction 
of quenchings from a fixed temperature that achieve the ground state varied considerably 
across system sizes. While the quenching of the energy of small systems occasionally resulted 
in achieving a ground-state structure, quenches of the large systems never resulted in a 
ground-state structure. The measure of the ground-state manifold becomes increasingly 
small as the system size increases. Given that the number of inherent structures increases 
exponentially with respect to system size,— it is unsurprising that larger systems would not 
sample the capture basins for ground states as frequently as smaller systems. System-size 
effects have been also been identified in the binary Lennard- Jones system. 25 Due to the 
system-size effects, we should not extrapolate true dynamical information from our results. 
Here, we present the analysis for 418-particle systems, noting that the trends across x values 
are similar to those in the 2340-particle system. 

The properties of the excited-state (or thermalized) structures behave as expected. The 
thermodynamic properties lagged the equilibrium values as the system was cooled, and for 
slower cooling rates, the lag in thermodynamic properties was smaller than that for fast 
cooling rates. However, we find that the rate at which the system was cooled provided no 
distinguishable effects on the properties of the inherent structures when compared across 
various cooling rates. This should be contrasted with conventional glass-forming systems 
such as binary Lennard- Jones where systems cooled under slow rates found deeper wells in 
the energy landscape than systems cooled under faster rates. 18 

In our previous research, we have identified the approximate range for transitions from 
crystalline or wavy- crystalline to liquids to be in the temperature range of T* < 0.00075. 2 
We find that the temperature at which the inherent structures were sampled from had no 
distinguishable effect on the thermodynamics properties or structures. As the system was 
cooled, the systems hopped across various capture basins without regard to temperature. 

Figure |3]exhibits the dynamics of systems as they sample inherent structures for x=0.6004 
and 0.8086 for 418 particles cooled at a rate of 7=10 -7 . The figure shows the potential energy 
of inherent structures sampled every 30 time steps. In the figure, there are two classes of 
behaviors. For \ = 0.6004, the system samples higher-energy inherent structures and also 
frequently samples ground states, those systems with e < 10~ 10 . For x = 0.8086, the system 
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samples mostly higher-energy structures and occasionally, though far less frequently than for 
smaller x, samples ground states. For both values of \i t ne system hops to a new inherent 
structure between samplings. There were no consecutive samples that were identical except 
for very low temperatures far below the melting point {e.g., T* = 0.00015). The rate at 
which the system samples capture basins does not appear to slow down significantly when 
the temperature is reduced as one might expect. This phenomenon occurred for all cooling 
rates we tested. In addition, even when sampling inherent structures every three time steps, 
the system never remained in the same basin for consecutive samplings. Note that the depths 
of the quenches for ground-state structures were deeper for \ — 0.8086 than for \ = 0.6004. 
This demonstrates the issue that arises with finite-precision computing because these scaled 
energies should vanish to zero in a mathematical sense. In a practical sense, systems with 
energies below this threshold and orders of magnitude below the clusters of systems higher 
up are ground states. 

We find that the structural features of the inherent structures did not vary with temper- 
ature either. Figure H] shows the structure factor for inherent structures with \ = 0.6004 
ensemble-averaged for several temperatures. However, when observing the structure factors 
for individual configurations, there is some small variability in the shape of S(k). How- 
ever, there are no clear correlations with temperature. In the figure, the structure factors 
associated with each temperature overlay nearly perfectly. Even at very low temperatures 
(T* = 0.00015), the structure factor for the inherent structures does not differ from those 
obtained from a liquid state. 

For familiar systems, as the temperature is reduced, the potential energy of the inherent 
structures typically becomes lower as the system continues to reduce the potential energy 
through structural rearrangements, sampling deeper and deeper basins. However, we find as 
the temperature is reduced, the energy of the inherent structures does not appear to decrease. 
This reveals that the energy landscape is relatively flat and insensitive to temperature. The 
heights of the barriers, which we explore in Sec. IVI1 are therefore expected to be relatively 
small. 

In Fig. [5], we plot the frequency in which the cooling and quenching procedure yielded a 
ground state for 418 particles. The plot shows the fraction of inherent structures that were 
ground states as a function of x f° r 3000 inherent structures for each \. We have found that 
the temperature from which the inherent structures are sampled does not have a significant 
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FIG. 3: Scaled potential energy e of inherent structures as a function of time for x=0.6004 and 
0.8086. The lower plots show the temperature schedules as a function of time. The temperature was 
reduced in a piecewise linear manner with an effective cooling rate of 7 = 10~ 7 . Inherent structures 
were sampled every 30 time steps. As a function of time, the systems continually hop from one 
capture basin to another since no two consecutively-sampled structures have identical energies. 
We found that the effective cooling rate did not have a significant effect on the fluctuations of the 
energies of inherent structures. At low temperatures, systems evidently are not trapped in deep 
energy wells. For x = 0.6004, a significant fraction of inherent structures have an energy below the 
ground-state threshold of e < 10 -10 , while only a few inherent structures lie below this threshold 
for x = 0.8086. 

effect on the rate at which the ground state is achieved. However, the value of \ affects 
the frequency at which the ground state is found. In general, as x is increased from 0.6 to 
0.9, the rate at which the ground states are found is reduced. However, for the "disordered- 
ground-state" region, x < 0.58, ground states were found infrequently. For example, in Fig. 
EJ the frequency in which ground states are found for x — 0.5622 is much smaller than that of 
X = 0.6004. This suggests that the relative availability of the ground-state manifold is small 
for x < 0.58, becomes larger in the wavy-crystalline region (x > 0.58), and then diminishes 
as the terminal x of 0.9 is reached. This shows that changes in x (and equivalently, density 
for fixed K), fundamentally affect the shape and size of the ground state manifolds. An 
expansion or compression of the system does not simply rescale the energy landscape while 
maintaining its overall shape. Systems containing 2340 particles never achieved ground state 
energies, presumably because of the exponential growth in the number of inherent structures 
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FIG. 4: (Color online) Structure factor for inherent structures with x=0-6004 and iV=418 and at 
various dimensionless temperatures T* - 0.00195 (circles), 0.00135 (squares), 0.0075 (diamonds), 
and 0.00015 (triangles). Each temperature is ensemble averaged over 50 configurations sampled 
from a molecular dynamics trajectory. Despite the range of temperatures, temperature has little 
effect on the inherent structures. 

with system size.— 

V. RESULTS: STRUCTURAL CHARACTERISTICS 

Figures [6] and [7J show inherent structures obtained via the cooling and quenching proce- 
dure. These structures display features that are common among inherent structures. Upon 
examining the inherent structures, we found several structural motifs - five-particle rings, 
wavy grain boundaries, and vacancy-interstitial defects. The frequency of these features 
varies depending on \- F° r 0.5 < \ < 0.6, the five-particle ring, as demonstrated in Fig. [6l 
is the most prevalent feature. The five-particle ring arises due to the local interactions of 
constituent particles. The nearest- and second-nearest neighbor distances lie within the first 
and second wells of v(r) respectively. The ratio of the relative distances creates a frustration 
that favors the local five-particle ring. Evidently, these local forces are more dominant than 
the long-range forces that promote stealthiness. 

For 0.58 < x < 0.89, a common feature of inherent structures is a wavy grain boundary 
as evident in Figs. |6] and [7J In these images, the grain boundaries manifest themselves in a 
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FIG. 5: Fraction of inherent structures obtained via the cooling and quenching procedure with 
7 = 1CP 8 and N = 418 that are ground states (i.e., scaled potential energy e < 1CP 10 ). Data 
represent the average of four independent cooling runs with 750 samples for each run. \ values of 
0.6004 and 0.7033 achieved the highest fraction of ground states. 

more wavy form than one might see in a Lennard- Jones inherent structure, in which particles 
prefer to align in straight lines. This x range is associated with wavy-crystalline ground 
states, which appear as a uniformly sheared triangular or square lattice. This behavior 
evidently manifests itself in the inherent structures where the grain boundaries appear to 
have some waviness. For these x values, the long-range nature of the potential is sufficiently 
strong to inhibit five-particle rings in favor of the wavy grains. Lastly, for x — 0.8995, the 
last image in Fig. [3, the inherent structures either have a large grain boundary of perfectly 
aligned particles, vacancy-interstitial pairs, or both. These inherent structures also exhibit 
polycrystallinity. In the cases of x > 0.58, the multiple wavy domains appear to be in 
structural conflict with each other. This is best demonstrated by Fig. [6] for x — 0.6004 
where several domains appear to meet in the middle of the system. This conflict evidently 
"jams" the system and prevents it from reaching the ground state. 

Observing the structure factor for inherent structures shows that inherent structures 
with smaller x are more stealthy and more hyperuniform than those at higher x- How- 
ever, inherent structures that are not ground states were neither completely stealthy nor 
hyperuniform. In Fig. [HI we show the ensemble-averaged S(k) for 50 inherent structures 
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FIG. 6: (Color online) Inherent structures produced from the cooling and quenching procedure 
sampled from a high-temperature liquid (T* = 0.0018) for \ = 0.5191 (left), 0.5622 (middle), and 
0.6004 (right) for 418 particles. These structures demonstrate the five-particle clusters that are 
common in inherent structures for these x values. The five-particle rings are highlighted by a 
surrounding circle. The buildup of wavy grains are visible for % = 0.6004. 
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FIG. 7: (Color online) Inherent structures produced from the cooling and quenching procedure 
sampled from a high-temperature liquid (T* = 0.0018) for % = 0.7033 (left), 0.8086 (middle), 
and 0.8995 (right) for 418 particles. These structures demonstrate the wavy grains and vacancy- 
interstitial defects common to inherent structures for x > 0.58. As x increases, the inherent 
structures are more ordered, though grain boundaries prevent the systems from minimizing the 
energy. 



sampled at T* = 0.00105 for various x- F° r the smallest x value, the small-k tail of the 
structure factor initially decreases and then rapidly increases as k is increased toward unity. 
The shape of the small-A; tail of S(k) remains constant but is shifted to higher values (less 
hyperuniform) as % increases to % — 0.5813. 
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However, as x leaves the disordered-ground-state regime, the shape of the curve changes, 
as exemplified by the curves for x=0.6004 and 0.7033. In Fig. [HI the curve for \ = 0.6004 
demonstrates the new shape. For the small-A; region, S(k) initially has a negative slope 
which quickly turns positive, and nearly flattens before diverging near K. The shape of this 
ensemble-averaged S(k) is due to the relatively large number of ground states that arise in 
the inherent structure analysis for this x- F° r X = 0.6004, nearly 40% of inherent structures 
are ground states, and this S(k) represents the average of the ground-state S(k) and the 
inherent structure S(k). A similarly large fraction of inherent structures are ground states 
for x=0.7033. If the ground states are removed from the ensemble-averaging, the structure 
factor appears more like those associated with x=0.8086 and 0.8995. For these \i S{k) has 
a initial negative slope that turns positive near fc=0.5, which then diverges near K. 

It is interesting to note that the way in which systems get "stuck" in an inherent structure 
for x m the disordered regime is different than that of the wavy-crystalline regime. For 
X < 0.58, the systems struggle most to minimize S(k) for k > 0.5. The length scale 
associated with k > 0.5 corresponds to local interactions in real space (i.e., on the order of 
a few particle diameters). This manifests itself with the five-particle rings. For x > 0.58, 
the systems have difficulty minimizing S(k) near zero and unity. The fc-values near zero 
represent longer-range interactions than the fc-values near unity. The inherent structures 
for x > 0.58 are therefore "stuck" on a longer length scale. It is no surprise then that the 
inherent structures appear locally ordered but are frustrated with a grain boundary. 

It is obvious that as x increases the extent of hyperuniformity of the inherent structures 
diminishes. This is shown in an ensemble sense in Fig. [HJ However, we have also shown 
this general relationship for individual structures. In Fig. [9j we display the hyperuniformity 
coefficient Co for the log-polynomial fit to S(k) in Eq. (fT3]) . In the figure, each data point 
represents one inherent structure. There are two clusters of points, where the cluster of 
points on the lower left are typically ground states where S(k) vanishes as k approaches 
zero, but due to numerical precision the hyperuniformity coefficient Co remains finite. In 
the cluster in the upper right of Fig. [9j there is a linear relation between Cq and log?]. 
Also, it is clear that as x increases, the hyperuniformity parameter Cq also increases. A few 
outlying structures are present for x=0.7003, 0.8086 and 0.8995. These structures appear as 
triangular lattices with very slight waves or shear applied to them, a sufficient amount that 
they are not ground states, but still have very small limiting values of S(k) as k approaches 
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FIG. 8: (Color online) Structure factor for inherent structures as a function of % ensemble averaged 
over 50 configurations sampled from a molecular dynamics trajectory at T* = 0.00105 for a system 
of 418 particles. As x increases, the inherent structures become less stealthy. Also, the extent of 
hyperuniformity decreases when x is increased. 

zero. 

The stealthiness of the inherent structures, as characterized by the metric rj, in general 
does not have a relation to the six-fold bond-order orientiational order of a system. Figure 
[T0l shows the bond-order parameter ^/q as a function of stealthiness 77 for various inherent 
structures. For each \ value, there is a range of available 77 and Initially, for \ near 
0.5, these ranges are rather narrow, and as x increases, these ranges become much broader. 
For large x=0-8086 there is a large diversity of and 77 available for inherent structures. 
However, for x=0.8995, there is a narrow band of available rj but a large range of ^6 
depending on the number and orientation of opposing crystalline domains. The group of 
structures for x=0.8995 and r\ ~ 0.02 are those structures with only vacancy-interstitial 
defects. 

The bond-order parameter also shows considerable variation when plotted against the 
hyperuniformity parameter Cq. Figure [TT] shows Co as a function of ^6 for various \ values. 
As x increases, the variation in grows while Co generally falls into a narrow range of 
values. This figure demonstrates that hyperuniformity does not necessarily require local 
six- fold ordering. Systems with x=0.7033 have the largest range of available local six- fold 
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FIG. 9: (Color online) Hyper uniformity parameter Cq versus stealthiness metric 77 for inherent 
structures generated from the MD and quench procedure. There is a clear relation between the 
hyperuniformity coefficient Co and the stealthiness metric rj, although there exist outlying config- 
urations. There is also a positive relation between Co and x- Those configurations in the lower left 
with r] < 10 -9 are typically ground states. 800 configurations are shown for each \ value. 



ordering for more hyperuniform structures because \l/6 ranges from as low as 0.1 to just 
below unity. 



VI. PATHS FROM INHERENT STRUCTURES TO GROUND STATES 

The nudged-elastic-band algorithm was used to determine the minimum-energy path 
from a ground-state structure to an inherent structure nearby in configurational space. We 
find that the height of the barrier between the ground state and the inherent structure 
generally increases with \. However, there are instances that deviate from this general 
situation. Figure [12] illustrates the minimum-energy paths found for various x from the 
ground state (point zero on the reaction pathway) to the inherent structure (last point on 
the reaction pathway). The pathway is discretized over fifty images. In all the cases, the 
landscape is locally quadratic near the ground state and inherent structure. The maximum 
for each pathway represents the saddle point in the energy landscape separating the ground 
state and inherent structures. Connections from the triangular lattice to inherent structures 
were challenging to construct because of the zero-energy modes present in triangular-lattice 
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FIG. 10: (Color online) Bond-order parameter versus stealthiness metric 77 for inherent struc- 
tures generated from the MD and quench procedure. As % increases, there is an increasing degree 
of six-fold ordering. However, for large Xi the range of available ^6 extends larger due to the 
polycrystalline nature of the configurations. 800 configurations are shown for each x value. For 
clarity, we omit those structures with rj < 10~ 5 so that the plot can provide more detail about the 
inherent structures. 

systems. The nudged-elastic-band algorithm encountered difficulty in traversing through 
zero-energy valley and then uphill toward a saddle point. In all cases except x=0-8995, we 
report transitions from non-lattice ground states to inherent structures. 

For x < 0.6 (left in Fig. fT2l) . the barrier height, e sa d is significantly large compared to the 
difference in the energies of the ground-state and inherent structures {i.e., e in h~e gs <C £ S ad)- 
When x becomes larger, as in the case for x=0.7033 and higher, middle and right in Fig. 
[T2| e sa( i is nearly equivalent in height. For the case of x=0.8995, the minimum-energy path 
appears to flatten out as it approaches the inherent structure. In reality, there is a maximum 
energy along the path, although it is difficult to discern due to the scaling in the plot. In 
the cases where the energy of the inherent structure is nearly equal to the barrier height, 
small perturbations from the inherent structure can allow the system to relax to the ground 
state. On the other hand, for smaller x, relatively higher-energy perturbations are necessary 
for the inherent structure to climb the energy barrier and fall to the ground state. While 
there are exceptions to this phenomenon, they usually occur when the inherent structure 
and ground state are not nearby in the energy landscape {i.e., there is a another inherent 
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FIG. 11: (Color online) Hyper uniformity parameter Cq versus bond-order parameter ^6 f° r inher- 
ent structures generated from the MD and quench procedure. Despite the broad range of six-fold 
ordering at larger Xi the level of hyper uniformity remains consistent across x values. 800 configu- 
rations are shown for each x value. For clarity, we omit those structures with Co < —12 so that 
the plot can provide more detail on the inherent structures. 
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FIG. 12: (Color online) Representative minimum-energy paths from ground states (point zero 
on the reaction coordinate) to inherent structures (last point on the reaction coordinate) along a 
generalized reaction pathway for N = 418 for (a) 0.5 < x < 0-6, (b) x = 0.6004 and 0.7033, and (c) 
X = 0.8086 and 0.8995. The local maximum corresponds to a saddle point in the energy landscape. 



structure closer to the ground state). Along these paths, the configurational proximity 
metric p monotonically decreases along the path from the ground state to the inherent 
structure. The metric p/N for the proximity per particle between the inherent structure 
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FIG. 13: (Color online) Configurations (top row) and associated Voronoi tessellations (bottom row) 
for the minimum-energy path from ground state (left column) through the saddle point (middle 
column) to the inherent structure (right column) for x = 0.5383. The transitions are from the set 
described in Fig. [I2j 

and the ground state is of the order 0.1. On average, particles are collectively displaced by 
10% of the nearest neighbor spacing to achieve a ground state. 

In Figs. H3]and[TH we display the ground state, saddle point, and inherent structure and 
the associated Voronoi diagrams for x=0.5383 and 0.7033. These are the same paths shown 
in Fig. [12j These systems are representative of the types of local rearrangements found in the 
transition from ground states to inherent structures. In Fig. [131 the ground state has wavy- 
crystalline characteristics where there is alignment in one direction. The Voronoi diagram 
appears to have highly ordered arrangement of polygons. In the saddle point images, one can 
see the appearance of the five-particle rings. In the inherent structure, the defect appears 
to be localized, however, the defects are more easily discerned in the Voronoi diagram. 

For x=0.7033, the buildup of the grain boundary is clear in Fig. [TH The ground state 
is wavy-crystalline and the saddle point shows the appearance of a vertical and horizontal 
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FIG. 14: (Color online) Configurations (top row) and associated Voronoi tessellations (bottom row) 
for the minimum-energy path from ground state (left column) through the saddle point (middle 
column) to the inherent structure (right column) for \ = 0.7033. The transitions are from the set 
described in Fig. [12j 



grain. The inherent structure shows the apparent mismatch between lines of particles that 
is responsible for the inability to reach the ground state. The Voronoi diagrams in this case 
visually display the source of the frustration in the system. While far from the grain, the 
polygons in the Voronoi digram for the inherent structure are quite regular and ordered, 
those near the grain boundary have much more variation in shape. For x=0.8995, the rear- 
rangements found in the transition from a perfect lattice to one with a vacancy-interstitial 
pair are very local (albeit collective), only involving the few particles immediately surround 
the interstitial- vacancy pair. In all cases observed, the rearrangements from ground states to 
inherent structures involve only local, but collective, rearrangements (within a few particle 
diameters). As highlighted here, they typically involve the formation of a grain boundary 
or a five-particle ring. 



27 



VII. PARTICLE REARRANGEMENTS FROM RSA PATTERNS TO GROUND 
STATES 



Configurations generated by random sequential addition near the saturation densities are 
known to suppress long- wavelength scattering, described as "ultratransparent,"— though 
images of the structure factor for the RSA configuration demonstrates that it is not hyperuni- 
form and not absolutely stealthy— From Ref. [34], as k approaches zero, the structure factor 
approaches a value of 0.059, corresponding to a value of the hyperuniformity parameter Co 
of -1.23. 

In Figs. [T5] and [TBI we display the initial RSA configurations and the final ground-state 
configurations as obtained from the collective coordinate approach, for x values of 0.0991 
and 0.3657. Each figure also displays the corresponding Voronoi tessellation derived from the 
point patterns. In both cases when observing the dynamics of the minimization algorithm, 
one observes a general spreading of particles due to the repulsion for small r. These collective 
rearrangements are very local and there are no large global rearrangements involved. While 
every particle is displaced collectively from their initial location, each particle is displaced 
by a small fraction of the mean-nearest neighbor distance, as measured by the proximity 
metric per particle, p/N. 

In Fig. [151 one can observe similar arrangements of particles in both the RSA configu- 
ration and the ground-state configuration. These images would be difficult to distinguish 
visually. Several structural features remain in tact aside from the general spreading of par- 
ticles. In the Voronoi tessellations of the RSA configuration, there are areas where the local 
density appears higher than in other areas. However, in the tessellation of the ground state, 
these areas of higher local density have been relaxed away. The density appears much more 
uniform in the ground state than the RSA configuration. However, the visual difference 
between the two images are subtle. In Fig. [161 the spreading remains local, although it is 
clear that the repulsion is stronger. The particles, on average, travel a larger distance and 
fewer structural features of the RSA pattern are maintained. The Voronoi tessellation of the 
ground state appears to have more uniform cell sizes than the RSA configuration, though 
the differences are subtle. 
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FIG. 15: (Color online) (Top) Configurations with p = 0.1953 and x = 0.0991 and (bottom) 
the associated Voronoi diagrams. The left images are the initial RSA configurations and the right 
images are the ground states. The trajectory from the RSA configuration to the ground state 
appears to arise from a gentle repulsion. The diameters of the particles correspond roughly to the 
assigned RSA diameter. The dark lines represent the system box. 
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FIG. 16: (Color online) (Top) Configurations with p = 0.05434 and x = 0.3657 and (bottom) 
the associated Voronoi diagrams. The left images are the initial RSA configurations and the right 
images are the ground states. The trajectory from the RSA configuration to the ground state 
appears to arise from a gentle repulsion. The diameters of the particles correspond roughly to the 
assigned RSA diameter. The dark lines represent the system box. 
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TABLE I: Comparison of RSA configurations to their associated ground states. 
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1.866 


5.928525xl0 4 


0.06614 


4.653 


0.00621 


0.19653 


0.0991 


737 


1.883 


5.462281 xlO 4 


0.05879 


5.258 


0.00713 


0.14660 


0.1364 


733 


2.155 


5.962653xl0 4 


0.07407 


4.915 


0.00671 


0.09893 


0.1995 


742 


2.639 


6.560097xl0 4 


0.07569 


6.354 


0.00856 


0.06851 


0.2902 


734 


3.154 


7.001039 xlO 4 


0.08924 


6.535 


0.00890 


0.05434 


0.3657 


741 


3.559 


7.501234xl0 4 


0.09812 


6.943 


0.00937 


0.04678 
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0.11910 


8.342 
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We have quantified the differences between the RSA structures and the ground states for 
several RSA configurations at various \ values. Table |T] displays the parameters of these 
systems as well as the differences in potential energy between the RSA system and ground 
state, the stealthiness metric for RSA configurations r] rsa , and the configurational proximity 
metric p. The stealthiness metric r] for ground states vanishes identically to zero. 

There are some variations in energies and metrics that are attributed to finite system 
size. Structural features of the RSA patterns can vary across samples, and these variations 
show up in Table [J particularly for small x- 111 general, the difference between the potential 
energy of the RSA system and the ground state increases uniformly with x, but is on the 
order of 10 -4 . The stealthiness metric also increases uniformly with \ as is expected, since 
decreasing the density of RSA configurations will rescale the structure factor to be less 
stealthy. The configurational proximity metric p increases uniformly with x, ranging on a 
per particle basis of 0.00621 to 0.01478 for these \- These metrics show that particles in 
the RSA configuration need to be collectively displaced locally by a small fraction of the 
mean-nearest-neighbor distance in order to produce stealthy, hyperuniform systems. 

VIII. DISCUSSION 

In this paper, we have shown that the energy landscape associated with the /c-space 
overlap potential is relatively flat and devoid of deep energy wells. This is to be contrasted 
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with Lennard- Jones systems, which possess rugged energy landscapes. The sampling of 
inherent structures from the energy landscape is independent of the temperature from which 
the system was sampled and the rate at which the system was cooled. Five-particle rings, 
which are related to the local interactions between particles, disrupt the ability of a system 
to become a stealthy ground state, while grain boundaries, which arise for \ values in the 
wavy-crystalline regime, are attributed to the longer-range interactions in a system. While 
the hyper uniformity parameter Co and the stealthiness parameter 77 are positively correlated 
with each other, the bond-order parameter \l/ 6 can display a broad range of characteristics 
for various x values. We have used the nudged-elastic-band algorithm to show that local 
particle rearrangements can allow systems to be perturbed over a relatively small energy 
barrier from an inherent structure to a ground state. These rearrangements correspond to 
a low configurational proximity metric p. Lastly, we have shown that highly local collective 
particle displacements, as quantified via the configurational proximity metric p, are sufficient 
to convert a RSA configuration into a stealthy, hyperuniform ground state. 

While these studies represent another contribution to understanding the nature of these 
unusual long-ranged potentials, there are some natural extensions to this work that can fur- 
ther illuminate the relationship between particle interactions, stealthiness, and hyperunifor- 
mity. We have used simple metrics Co and 77 as measures of stealthiness and hyperuniformity, 
there are a multitude of other metrics that can be constructed. Measuring hyperuniformity 
for a single finite system is challenging due to the fluctuations and noise that can arise in the 
small- k behavior of S(k). Ideally, one would prefer to have a large system where the number 
variance in an observation window can be measured accurately. However, we are limited 
in system size because of numerical methods. While our measurement of stealthiness is 
directly related to the structure factor, we assumed equal weights for all the S(k), while one 
could assign weights differently to favor small-A; scattering or near- if scattering. Currently, 
there are no obvious metrics that stand out as superior to those used here. However, the 
development of new, improved metrics remains a potential direction of future research. 

While our paths connect ground states to inherent structures, a more generic algorithm 
taking inherent structures to ground states would be valuable. For instance, our methods 
relied on knowledge of the ground state and inherent structure a priori. However, a more 
useful method would be to find the nearest ground-state structure given only an inherent 
structure. This algorithm would be far more powerful than the current method. One could 
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further quantify the possible rearrangements in the inherent structures in order to help 
produce such an algorithm that relies on trial displacements that search for small energy 
barriers nearby. This could lead to a statistical analysis of the types of rearrangements made 
in a system to determine which are more likely to lead to a ground state. In addition, a 
more in-depth study of the rearrangements from RSA configurations to hyperuniform and 
stealthy configurations would be meaningful. Perhaps one can distinguish the relative con- 
tributions of the features of the interaction potential that contribute to the hyperuniformity 
and stealthiness of a system. 

The precise role of the pair interactions in the formation of disordered ground states and 
stealthiness is unknown. Here, we have determined that five-particle rings are a consequence 
of the first two minima in the overlap potential v(r). The ground-state behavior holds for a 
class of potentials 3 and we strongly suspect that similar behavior in the inherent structure 
analysis may generally hold as well. A systematic study involving the determination of 
ground states for the truncated overlap potential in Eq. (J^J) might help understand how 
many local minima in v(r) are necessary for disordered ground states to arise. The potential 
could be truncated at the first, second, third, etc. minimum and numerical procedures could 
be used to find ground-state structures. This would encounter issues of rigor in that the 
structures found by numerical methods may not touch the lower bound on 0. The trivial 
lower bound for <fi is zero (due to the nonnegativity of the overlap potential), but for the 
densities associated with the disordered and wavy-crystalline regimes, the improved lower 
bound, that being S(k) must vanish for all k < K, is not applicable to the truncated /c-space 
overlap potential. 

In addition, we observed that the shape of V(k), which assigns weights to S(k) in the 
potential-energy function, apparently influences the shape of the function S(k) for the in- 
herent structures for x < 0.58. Above this x value, there was no obvious connection between 
V(k) and S(k) for the inherent structures. This raises the question as to whether the shape 
of V(k) actually influences the shape of the function S(k) for inherent structures. A sim- 
ple study varying the shape of V(k), but maintaining the compact support at K and the 
positivity of V(k) would help answer this question. 

Lastly, while this collective coordinate setup was used understand and design new ma- 
terials in two-dimensions, such as photonic bandgap materials,— the extension to three 
dimensions should be particularly fruitful. While we have some understanding of stealthy 
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ground states in three dimensions,- 1 ^ the relations between inherent structures and ground 
states may be different in three dimensions than two dimensions. Connecting these struc- 
tural rearrangements and physical properties can help to contribute to the next class of 
three-dimensional photonic bandgap materials. 
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